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We investigate the two-magnon Raman scattering from an anisotropic S =1/2 triangular Heisen- 
berg antiferromagnet CS2CUCI4. We find that the Raman response is very sensitive to magnon- 
magnon interactions and to scattering geometries, a feature that is in remarkable contrast with the 
polarization-independent Raman signal from the isotropic triangular Heisenberg antiferromagnet. 
Since a spin-liquid ground state gives rise to a similar rotationally invariant Raman response, our 
results on the polarization dependence of the scattering spectrum suggest that Raman spectroscopy 
provides a useful probe, complementary to neutron scattering, of the ground-state properties of 
CS2CUCI4, particularly whether the time-reversal symmetry is broken in the ground state. 



I. INTRODUCTION 

Recently Heisenberg antiferromagnets on the trian- 
gular lattice have attracted considerable experimental 
and theoretical interest. Among them, CS2CUCI4 has 
been under particular scrutiny as it provides an inter- 
esting example of a spatially anisotropic spin-1/2 trian- 
gular antiferromagnet Much of the interest in this 
compound stems from its unusual, non-classical mag- 
netic properties, arising from the competition between 
the spatial anisotropy, Dzyaloshinskii-Moriya (DM) in- 
teractions, and quantum fluctuations. 

Extensive neutron scattering studies^ on the mag- 
netic properties of CS2CUCI4 revealed several inter- 
esting features. First, despite frustration and low- 
dimensionality, a long-range magnetic order develops at 
low temperatures: the observed spin order is incommen- 
surate and sets in at temperatures below TV = 0.62 K. 
Magnetic excitations above this ground state are also 
quite unusual. While the low-energy excitation spectrum 
contains well-defined sharp modes, as expected for an or- 
dered state, a broad continuum is formed at intermedi- 
ate and high energies. A number of theoretical proposals 
have been made to explain the origin of this continuum. 
It has been suggested that the existence of a continuum is 
an indication that the system is proximate to a spin liquid 
phase that determines the behavior of excitations except 
for low energies^ - — £ An alternative suggestion is that the 
continuum might originate from magnon-magnon scat- 
tering which is enhanced in non-colinear magnetsi&i 

In view of the ambiguity in the interpretation of the 
neutron scattering data, a complementary experimen- 
tal analysis of the magnetic properties of CS2CUCI4 by 
a different technique is highly desirable. A very ef- 
fective and frequently used experimental tool to study 
low-temperature properties of low-dimensional quantum 
magnets is the two-magnon Raman scattering.— ~— The 
two-magnon Raman intensity is directly related to the 
spectrum of two interacting magnons in a total spin zero 



state at vanishingly small momentum and weighted by a 
form factor that is dependent on the polarization of the 
incident light. It contains detailed information on the 
two-magnon density of states and the magnon-magnon 
interactions. Therefore, direct comparison of experimen- 
tal spectra with those obtained from theoretical anal- 
ysis can lead to rather accurate estimates on superex- 
changc and DM interactions. In addition, the analy- 
sis of the polarization dependence of the magnetic Ra- 
man scatterin g 16 ! 17 might shed some light on whether 
the ground state is ordered, as neutron scattering exper- 
iments have suggested, or is actually in a spin liquid state, 
as some theories suggest. A pronounced polarization de- 
pendence would indicate that magnetically ordered state 
is the most probable candidate for the ground state, and 
that the observed continuum is due to relatively strong 
interactions between magnons. On the other hand, if 
the Raman scattering depends weakly on the scatter- 
ing geometry, the continuum in neutron scattering might 
be due to unconventional excitations above a spin liquid 
ground state. 

In this paper, we carry out a theoretical analysis 
of two-magnon Raman scattering from an anisotropic 
S = 1/2 Heisenberg antiferromagnet (HAF) on a 
triangular lattice (see FigfT]). To evaluate the two- 
magnon Raman spectra, we use the well-established, 
semi-phenomenological Loudon-Fleury (LF) approach^ 
By chosing model parameters relevant for the com- 
pound CS2CUCI4, we find that the spectral shape as a 
function of frequency is sensitive to 1/S corrections of 
the magnon spectrum and is strongly modified by the 
magnon-magnon interactions in the final state. The in- 
tensity of the two-magnon peak also significantly depends 
on the scattering geometry which is in contrast with the 
polarization independence of the magnetic Raman re- 
sponse in isotropic HAF on triangular latticed 

The paper is organized as follows. In Sec. [Hi we present 
a model of CS2CUCI4 and discuss its classical ground 
state, which is an incommensurate spin spiral with a 
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FIG. 1: Triangular lattice with anisotropic spin exchanges. 



pitch vector determined by the competition of anisotropic 
nearest-neighbor interactions. In Sec. IIII1 we first review 
results of the one magnon excitations in the anisotropic 
5 = 1/2 HAF to first order in 1/S. We show that al- 
though 1 / S'-corrections are present in the whole Brillouin 
zone (BZ), they are less drastic than in the isotropic case 
of the triangular lattice. This is due DM interactions, 
which suppress quantum fluctuations and open a gap at 
the ordering vector. In Sec. IIV1 we first review the LF 
formalism and then use it to calculate the Raman spec- 
tra at various levels of approximation, i.e., using only 
the bare magnon dispersion, using a magnon dispersion 
renormalized to order 1/S, and with final state interac- 
tions included. We show that the Raman profile is very 
sensitive to the magnon-magnon interactions and to the 
scattering geometry. Finally, Sec. [V] presents a summary 
of the work. 



II. MODEL 

We start with the following spin-1/2 model Hamilto- 
nian for Cs 2 CuCi4: 



H 



Si ■ Sj + Djj • (Si x Sj)] . 



(1) 



where (ij) refers to nearest-neighbor (NN) bonds on the 
triangular lattice, and and Dy are the symmetric and 
antisymmetric exchange constants. The antisymmetric 
spin exchange originating from the relativistic spin-orbit 
interaction is also known as Dzyaloshinskii-Moriya (DM) 
interaction. For CS2CUCI4, it is customary to denote the 
exchange constants along the horizontal bonds, which 
form quasi-one-dimensional chains, as J, and along 
the zigzag bonds as J'.— £ The DM interaction vanishes 
along the quasi-lD chains, whereas on the zigzag bonds, 
the DM vectors arc perpendicular to the triangular plane 
T>ij = ±(0,0, D) (Fig. [T|. Experimental measurements 



in high magnetic field have resulted in J w 0.374 meV, 
J' w 0.128 meV, and D = 0.02 meV^ 

The classical ground state of Hamiltonian (TTJ) is given 
by a spin spiral Si/S = cos(Q • r,) b + sin(Q ■ r^) c, where 

the pitch vector Q = Qh depends on the ratio J' /J. 
In the isotropic case J' = J, the ground state is the 
well known 120° non-collinear magnetic order with Q = 
27r/3i 19 ' 20 To understand the classical ground state in 
detail, we consider the energy of the magnetic spiral: 
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Minimization with respect to Q leads to the following 
three types of long-range magnetic order: (1) At J' > 2 J, 
the magnetic spiral reduces to collinear Neel order with 
ferromagnetic ordering of the spins along the chains. 
Along the c direction, spins on adjacent chains are an- 
tiparallel to each other. (2) At < J' < 2 J, the pitch 
of the spiral is given by Q = 2 arccos(— jj), which varies 
from 2ir — > tt as we vary J' from 2 J to 0. (3) For J' = 0, 
the system degenerates into decoupled antiferromagnetic 
chains. 



III. LARGE-5 EXPANSION 

The large-S* expansion about the classical spiral or- 
der can be significantly simplified with a locally rotated 
frame of referenced The spin components Si in a lab- 
oratory frame are related to those in the rotated local 
frame through 



Sf 

sv 

Sf 



S? cos Q — Sf sin Q, 
Sf sin Q + Sf cosQ 



(5) 



The spiral viewed from the rotated local frame corre- 
sponds to a simple ferromagnetic order Si = Si. We 
then employ the Holstcin-Primakoff transformation: 



st 

s- 



(0) 



(28-ata t ) 1/2 a, 
at(2S-afa i ) 1 / 2 



The magnon operators a\ and describe excitations 
around the spiral ground state. As we intend to study 
magnon interactions to first order in 1/S, we need to ex- 
pand the Hamiltonian in Eq. (TTJ) up to quartic order in 
the boson operators: 



H = E Q + 3JS(H 2 + H 3 + H 4 ). 



(7) 
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Introducing Fourier transform a, = X)k a k e /vN, 
the explicit expression of the various terms in the Hamil- 
tonian reads 



Ho 



ak a. 



-k + a], a f _ k ) 



(8) 



H 3 = l\l^sT,( C i+C 2 ) (4axa 2 - 44a 3 ), (9) 
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Here 1 • • • 4 denote ki • • • k 4 , and the summation in H3 
and i/4 is subject to momentum conservation module a 
reciprocal lattice vector: J^i k» = mod G. The follow- 
ing functions are introduced: 



2 ("^Q+k 



^Q-k) 



9 7 T 



5k — 2 (^Q+k "T" "'Q-kJ 

Here, for convenience, we rescale the interactions Jk and 
J k with respect to J, which is assumed to be J = 1. 

The quadratic Hamiltonian H 2 can then be diagonal- 
ized by a Bogoliubov transformation: 



r — j t i T 

^ k — J Q+k "'Q-k- 



(11) 



w k c k 
«kC k 



u k cl k 

«kC_ k 



(12) 



where c k ^ are operators for Bogoliubov quasiparticles. 
The coherence coefficients 



u k 



^k + £k 

2£ k 



_Bk 

\B k \ 



A k -E k 
2£ k 



satisfy u k — u k = 1, and 



Eu = 



A 2 
A k 



(13) 



(14) 



describes the quasiparticle dispersion (the energy of the 
quasiparticles is SJSE^). The diagonalized Hamiltonian 
H 2 is given by 



where 



H 2 = E 2 (Q) + Y,E k cic k 



£ 2 (Q) = ^(Ak^ + BkWkVk) 



(15) 



= -NJ A r 
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FIG. 2: Top: Renormalized magnon dispersion. Dotted line 
corresponds to the linear spin- wave dispersion E^, while red 
and blue solid lines correspond to the real and imaginary part 
Re(Im)Eu,, respectively, computed on a lattice of 252 x 252 k- 
points with artificial line broadening of r\ = 0.003. Middle and 
Bottom: 3D-plot of the ReE^ and of the ImE^, respectively. 
The spectrum is computed for J = 1, while other parame- 
ters describing interactions is CS2CUCI4, are correspondingly 
rescaled. All energies are measured in units of 3JS. 



gives 1/5 correction to the classical ground state energy 
Eq(Q). The 1/S correction to the ordering wave vector 
Q is determined by minimizing the sum 



E + E 2 = 3JS(S - l)NJl + 



3JS 



(17) 
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with respect to Q. The quantum correction is given by where E^ 4 '(k) is self-energy correction of order 1/S: 



AQ 



1__ Ly 



d 2 J 1 /dQ 2 N 



Ak — i?k ^Q+k 
2£ k dQ~ 



, (18) 



where Qo = 2arccos(— ^j) is the pitch of the classical 
ground state. 

The 1/S contribution from the quartic Hamiltonian 
can be obtained through a mean-field decoupling of H4. 
We first define 

G(k) = (4o k ), F(k) = (a k a-k> = (a k a -k>- ( 19 ) 

The quadratic Hamiltonian plus the decoupled H4 can 
be expressed as: 



A k 4 a k + — (a k a„ k + 4 a -k) 



where 



(20) 



= Av 



k-q + ^k-q 



Ak - Aq ) G(q) - ( B q + ) F(q) ,(21) 



q 

+ (^ k _ q + S k _ q -^-^)F(q) 



The magnon spectrum rcnormalized by the quartic 
Hamiltonian H4 becomes 



JAI -Bl = E k + E k 4) + 0(1/ S 2 ) , (23) 



,(4) 



k - 4 4 Vk+^)+25r« 

= (A k A k 4) -B k B^)/E k . (24) 

To obtain the 1/S correction from the cubic Hamil- 
tonian H3, we follow Ref. and consider interactions 
between quasiparticles c, : 



1/3 = i V2ibE[ $ i( k i' k2 ; k 3)44 2 ck 3 



-^2(ki,k 2 ,k 3 )c ki 4 2 c k 



h.c. 



(25) 



The vertex functions are given by (for simplicity, we de- 
note 1 = ki, 2 = k2, etc.) 

•,(!,*»)- * 2 (l,2,3) = l4y#.(26) 

where 
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$i(l, 2; 3) = Gx/^C/^/f + /- 7i J ) 
+G 2 /i 2) (/f/| 1) + /i 3) / a) ) 
+G 3 /i 3) (/| 1) /| 2) -/ W /i 2) ) 



$2(1,2,3) = Ci/^C/f /| 3) - /i 2) £ 3) ) 

+ C2/ ( 2 ) (/ (3) / (l)_ / (3) / ( 1)) 

+G 3 /i 3) (/! 1) /| 2) -/i 1) /i 2) ) 

and /i a) = VA* ± B a for a = 1, 2, 3. 

The triplic contribution to the self-energy is 



(27) 
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ki+k 2 = -k 



?ki + E^ 2 + E k +ir] 



(29) 



The first term in Eq. (|29| describes a virtual decay of 
a magnon into two-particle intermediate states. Terms 
with three creation(annihilation) operators, as will be- 
come clear in the next section, play no role in evaluating 
the magnon interactions within the Raman response. We 

(3) 

also note that S k is computed within so-called on-shell 
approximation. In this approximation the self-energy is 
evaluated at the bare magnon energy E k . 

Finally, the magnon energy renormalized by both the 
quartic and the triplic terms is given by 



Eu 



,(4) 



.(3) 



0(1/S 2 



(30) 



In Fig. [5] we plot the renormalized magnon spectrum 
i?k for parameters relevant to CS2CUCI4. One can see 
that the renormalization of the spectrum is much less 
pronounced than in the case of the isotropic triangular 
lattice ^ 1 ' 22 Moreover, the imaginary part of the magnon 
energy, ImE k almost vanishes in the whole BZ. Thus, the 
life time of the quasi-particles is very large at almost any 
momentum. This happens because the DM interaction 
opens a gap at the ordering Q vector, which significantly 
suppresses quantum fluctuations^ We note that, as cor- 
rections to the ordering vector Q, determined by Eq. (|18|) . 
are small at those values of D relevant to our calculation, 



we use the bare value of the ordering wave vector Qo for 
the remainder of the paper. 



IV. RAMAN INTENSITY 

A. Loudon-Fleury formalism 

Here we present the analysis of the two-magnon Raman 
scattering from the anisotropic triangular lattice. We 
employ the LF approach which models the interaction of 
light with spin degrees of freedom. The LF scattering 
operator is given by the photon-induced super-exchange 



operator* 



R 



18,23 



(£in ' <^i)(£out 



(31) 



weere 5^ denote basic vectors of triangular lattice: S± = 

(1,0), S 2 = (i,^) and <5 3 = (-|, %. J M defines 
the Heiscnberg exchange on the bond 5^. Since J M 
is anisotropic, the C^ v symmetry of the triangular lat- 
tice is broken. Thus, instead of using C3i,-irrcduciblc 
representations (A\, A 2 and E) for characterization of 
the polarizations, we determine polarizations of incom- 
ing and outgoing light as £j n = cos6>x + sin^y and 
£out = cos</>x + sinewy, where 9 and <f> are defined with 
respect to the x-axis. 

In terms of Bogoliubov quasi-particle c-operators, the 
LF scattering operator (|3"Tj) takes the following form: 



R = Y. M k( c kC-k + 4 c -k) = r ~ 

k 



(32) 



where Mk is bare Raman vertex, which is determined by 
the the magnon spectrum and by scattering geometry. 
The expression for Mk is given by 

M k = F 1 (k,e,<t>)u k v k + F 2 {k,e,(f>)(ul + vl) , (33) 

where we introduced the following notations: 



Fi(k,M) = 2SEiU/M(M)f,* 



(34) 



I(£i) 0.6 




FIG. 3: Bare Raman intensity (solid lines) and Raman inten- 
sity including the one-magnon renormalization of the spec- 
trum (dashed lines) at different polarizations of light de- 
scribed by <j>, and 9 scattering angles. Number of k-points: 
252 x 252. The broadening parameter is r\ — 0.003. is 
measured in units of 3 JS. 



The functions /^(fl, 0) = (ejjj-5 M )(e ut are symmetry 
weighting factors along the three basic vectors 6^ of the 
triangular lattice. 

The Raman intensity is obtained from the Fermi's 
golden rule, which on the imaginary frequency axis reads 



Jo 



dr < 



(T t (R(t)R))] 



(36) 



By analytically continuing Matsubara frequencies uj m = 
2'nmT onto the real axis as iu) m — > Q + ir], we obtain 
the Raman intensity as a function of the inelastic energy 
transfer = u>- m — w ou t of the incident photons. For the 
rest of the paper we assume the temperature T = 1//3 to 
be zero. 

To order 1/S, the Raman polarization operator 
(T t (R(t)R)) contains only terms like (T T (r+(T)r-)} + 
(T T (r~ (r)r + )) , where ^ are specified in Eq. ([32]) and 
Fig. EH). By Hcrmitian conjugation, it is sufficient to 
calculate J(r) = (T T (r~ (r)r + )) , which is depicted in 
Fig. [5b). 
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(35) 



B. Raman intensity without final state interactions 

We now focus on the contribution from two-magnon 
Raman scattering, which can be computed at different 
levels of approximation. We begin by first using only 
the bare spin-wave dispersion. Then we continue by in- 
cluding renormalizations of the one-magnon spectrum to 
1/S order. Calculating the Raman polarization operator 
(T t (R(t)R)) with the bare propagators of the Bogoliubov 
quasiparticles, G(k, iw n ) = l/(iu) n — E^), where E^ is the 
bare magnon spectrum gives the following expression for 
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FIG. 4: Variation of the maximum of the bare and renormal- 
ized intensities as a function of the scattering angle 9 com- 
puted at = 0.8 and Q = 0.6, respectively. 



the Raman intensity: 



7(Q) 



d 2 k / duoMl 
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1 



= — 277 9x1 
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k n - 2E U + ir) 



(37) 
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Fig. [3] shows the bare Raman spectra (solid lines) as 
functions of the transferred photon frequencies Q for 
three scattering geometries: (i) 9 = tt/2, = 0, (ii) 
9 = 2tt/3, = and (hi) 9 = 0, = 0. For two of the 
polarizations (i) and (ii) , the Raman spectra show a sim- 
ilar profile and intensity: the Raman response exhibits a 
peak at ~ 0.8 — 0.9, and the location of this peak corre- 
sponds to the twice the energy of the dominant van- Hove 
singularity of the one-magnon dispersion, cf. (Fig. [5]). 

The small Raman intensity at low energies is due 
to two reasons. First, the magnon energy is gapless only 
at the zone center and has a gap caused by the DM in- 
teraction at the ordering wave vector Qo- Thus, only 
magnons with momentum near k ~ can be excited by 
photons with frequency ft < 2E(Qq). As one can see 
from Fig. [2J the rather steep magnon spectrum in the 
vicinity of zone center gives rise to a small density of 
states. Second, the form of is such that it selects 
mostly wavevectors k ~ tt where the gap resides. 

In the # = = geometry, the Raman response is 
non- vanishing but very small for the whole energy range 
compared with intensities observed in other geometries. 
The relative smallness can be understood by compar- 
ing this result with the well-studied case of the square 
lattice, for which the LF operator in the A\ g geometry 
commutes with the Heisenberg Hamiltonian, and, as a 
result, the Raman response vanishes. In the case of the 
anisotropic triangular lattice, a non-commuting part of 
the LF operator (|31[) remains non-zero even for the A lg 
geometry. This part leads to small Raman intensities, 
scaling with the ratio (J 7 / J) 2 . Indeed, as we have dis- 
cussed earlier, the anisotropic triangular lattice can be 



viewed as an interpolation between the square and the 
isotropic triangular lattice by varying J'/ J. 

In Fig. SI the two-magnon peak is shown MS cl func- 
tion of the scattering angle 9. A strong dependence on 
the scattering geometry can be seen: the largest peak 
intensity is observed at the 9 = 7r/2,0 = 0, i.e. cross- 
polarization, which gradually decreases and reaches its 
minimum at 9 = or it (Though not exactly, the peak 
intensity is roughly proportional to sin 2 6*). This polar- 
ization dependence is consistent with the fact that the 
anisotropy in the present case resembles more that of a 
rectangular than of an isotropic triangular geometry. 

An angular- dependence analysis of the Raman spec- 
trum on CS2CUCI4 is highly desirable as it could pro- 
vide a potential diagnosis of whether a long-range spin 
order develops in the ground state. If an angular de- 
pendance similar to the one described above is observed, 
then the ground state is likely to be an incommensurate 
spiralj^i On the other hand, a Raman response that is 
independent of the scattering geometry might indicate 
a spin-liquid ground stated as proposed in some re- 
cent theories^— £ Here some precautions are necessary. 
Since the underlying Hamiltonian is spatially anisotropic, 
the spinon excitations of the proposed spin-liquid phase 
might inherit the crystal anisotropy to some degree, 
hence also giving rise to a polarization-dependent Ra- 
man signal. This question should be further studied but 
we believe that even if this is the case, the polarization 
dependence of the Raman response will be much weaker 
than in the case of spiral magnetic order. 

Next we incorporate 1/S corrections to the Raman 
spectrum. This can be easily done by replacing the bare 
energy with the rcnormalizcd magnon energy E^ (30) in 
the propagator G(k, iui n ) = l/(iuj n — E-k). The rcnormal- 
ized Raman spectra are shown by dashed lines in Fig. [3] 
We can see that in the renormalized spectrum the two- 
magnon peak appears at the energy Q ~ 0.6, which is 
slightly lower than the peak in the bare spectrum. This 
is in contrast with the case of isotropic triangular lattice, 
where the peak is shifted to higher energies. Once again 
it shows that the Raman spectra on the anisotropic tri- 
angular lattice has features of the both triangular and 
square lattices. 



C. Raman intensity with final state interactions 

Next we consider the final-state magnon-magnon in- 
teractions. Usually these are not small, particularly 
for S = 1/2. The effect of the final-state magnon- 
magnon interactions can be taken into account by 
computing vertex corrections to the bare Raman ver- 
tex. Here we consider only the leading 1/5 correc- 
tions. In this approximation, the vertex corrections can 
be obtained from an infinite summation of ladder dia- 
grams. These ladder diagrams are shown in Fig. ^jp) 
in terms of the two-particle (ir)reducible Raman vertex 
(7(k, p, cj n ,u; )) r(k, ui n , oj m ). These are related by the 



7 



Bcthc-Salpctcr equation: 

r(k,w„,w m ) = r~(k) + ^ 7(k, p,w„,cj ) 
G(p, w + w m )G(-p, -w )r(p, w , ui m ) . (38) 



The two-particle irreducible vertex can be decomposed as 
in Fig.[SJi): 7(k, p, w„, u ) = j 3 (k, p, u n , u a ) + 7 4 (k,p). 
The quartic vertex 74 (k, p) is identical to the two- 
particle-two-hole contribution from the H± term; its ex- 
plicit expression is given by 



74(k,p) = ^[4(A 



B + A p+k - -Bp+k + A p - Ak)MpU k Wpi'k + 



(/lp_ k - B p _ k + A p+k - B p+k - A p - A k ) (u 2 v ul + vpu£) - 
(2B p + B k ) (u 2 p + vl)u k v k ~ (2B k + B p ) (u k + vl)u p v f 



(39) 



where the functions A k and _B k are given by Eqs. (jlip . intermediate propagator, and can be written as 

The triplic vertex 73 (k, p, u>„, uj ) is obtained from the 
product of two vertices of the cubic term H3 and one 



73(k,p,w n ,a; ) 



1 

325 



[*i(p, k - p; k)$J (p, k p; k)G°(k - p, ilo 



$1 (-P, P k -k)<^ (-p, p k; -k)G°(p - k, iuj n - i0J o )] 
I 



(40) 



where the functions $i(l,2;3) and their complex conju- 
gates are given by Eqs. ([2^ ]) - ([2T| . To keep 73 (k, p, uj n , ui ) 
to leading order in 1/5, we retained only the zeroth order 
propagators G° for each intermediate line. We further 
simplify the expression (|40[) by assuming that the dom- 
inant contribution to the frequency summations in the 
Bethe-Salpeter equation ([3"81 comes from the mass-shell 
of the intermediate particle-particle propagators. This 
corresponds to the substitution of the intermediate fre- 



quencies by 



E p . The simplified 



expression for the triplic vertex then reads as 



73(k,p) 



1 

325 



$i(p, k - p; k) <f>{ (-p, p k; -k) 



2-Ek-p 



(E k -E p ) 2 -E£_ 



(41) 



The triplic vertex then depends only on momenta. 

Next we perform the frequency summation over u> on 
the right hand side of Eq. (|3"8")l as well as the analytic 
continuation iuj m — > D, + it] = z. With this, the reducible 
vertex T in the latter equation turns into a function of p 
and z only, leading to 



^L k , p (z)r p (z)=r-(k) 
p 

7(k,p) 



£k, P 00 = 5^ p 



z-2E„ 



(42) 
(43) 



which is an integral equation with respect to momentum 
only. Finally, the expression for the Raman intensity can 
be written as 



j(Q) ~ [j(n)-j(-n)] 



where 



J(Q) 



E 



irj 



2E U 



(44) 



(45) 



The fully renormalized intensity for the polarization 
with 9 = 27r/3, <f> = is shown in Fig. |6K. Despite the 
damping by the vertex corrections, the two-magnon peak 
survives and is further shifted towards lower energies. 
Apparently, the damping of the peak is less pronounced 
compared with the isotropic triangular lattice case.— At 
the higher energies, we also see the appearance of a broad 
continuum. We would like to point out that both the 
peak and the broad continuum are observed in geometries 
with almost perpendicular £i n and e ou t- Although the 
width and intensities of the peak and the profile of the 
continuum vary with the angle 9, the position of the two- 
magnon peak and the center of the continuum do not 
change significantly for 9 ~ tt/2. On the other hand, 
the Raman signal is extremely weak for nearly parallel 
geometries (9, <j> ~ 0), similar to the case without vertex 
corrections (see Fig. [3]). 

In order to disentangle the contributions coming from 
the triplic and quartic terms, we also compute the Raman 



R=> + < 



;(0) 



r h 
a) 




-k, -ico 



n b) 



c) 



co n +co m co +co m 
-ro n -<D n 




FIG. 5: a) Bare Raman vertex R from Eqn. (J32J) ; b) Raman 
susceptibility bubble; c) The integral equation for the dressed 
Raman vertex F in terms of the irreducible magnon particle- 
particle vertex 7; d) Leading order 1/S contributions to 7. 



spectrum with an irreducible vertex which includes only 
the quartic part. The comparison between the Raman 
spectrum computed with the full vertex and with the 
one containing only 74 (k, p) is presented in Fig. One 
observes that vertex corrections due to the quartic term 
split the sharp two-magnon peak into two peaks of com- 
parable intensities, which are, however, about the half 
of the intensity of the two magnon peak with final-state 
interactions. The triplic term modifies these two peaks 
quite differently. The lower energy peak is only weakly 
rcnormalizcd by the triplic term, while the higher energy 
peak is damped more strongly and is transformed into 
the broad continuum. 

Finally we note that the direct comparison of these re- 
sults with Fig. [3] should be taken with a certain precau- 
tion, since the artificial line broadening in Fig. [5] is larger 
by one order of magnitude. This is a consequence of a 
factor of 16 less k-points used in the latter case. This is 
because the kernel Lk, P (z) in the integral equation (|4"2")l 
and (H3J) is not sparse and has rank iV 2 x N 2 . Con- 
sequently, a moderate lattice size gives rise to a rather 
large dimension for the kernel. In the above calculations, 
we have chosen N = 69, leading to a 4761x4761 sys- 
tem which we have solved 100 times to account for 100 
frequencies in the interval f2 <E [0,2]. We also note that 
the kernel has points of singular behavior in (k, p)-space, 
which stem from the singularities of the Bogoliubov fac- 
tors and from the energy denominators in vertex func- 
tions. Here, we have chosen to regularize these points 



— Bare bubble 

— Renormalized bubble 

1 \ ■ ■ Renormalized bubble with vertex corrections 

\ 1 

\ 




9=2^/3, 6=0 



► Renormalizalion only by quartic term 

> Renonrjalizalion by both quartic and triplic terms 




FIG. 6: Effect of final state interactions on Raman inten- 
sity for 9 = 2ir/3,(f) — 0. Top: Bare intensity, an intensity 
computed with renormalized magnon energies and intensity 
computed with included final state interactions are shown by 
solid, dashed and dotted lines, respectively. Number of k- 
points: 69 x 69. The imaginary broadening is rj = 0.03. Q is 
measured in units of 3 JS. Bottom: Comparison of the intensi- 
ties computed with corrections only due to the quartic vertex, 
74 (k), and with corrections due to full vertex, 73 (k) + 74(k) . 



by cutting off eventual singularities in Lk,p- This can be 
justified because the weight of these points is negligibly 
small compared with the total number of points in the 
BZ. We have checked that this regularization does not 
significantly effect the obtained spectra. 



V. SUMMARY 

In summary, we have studied the two-magnon Raman 
scattering in the anisotropic triangular Hciscnbcrg anti- 
fcrromagnct considering various levels of approximation 
within a controlled l/S'-expansion. We have shown that 
the Raman profile is sensitive to the magnon-magnon in- 
teractions and to the scattering geometry. The calcula- 
tions indicate that the main effect of the magnon-magnon 
interactions is on the shifting of the two-magnon peak to- 
wards lower energies and on the formation of the broad 
continuum at the higher energies. We have also shown 



9 



that through exchange and DM interactions, the spatial 
anisotropy of the lattice is transferred to the magnon dis- 
persion of the spiral magnetic order. This makes the two- 
magnon Raman scattering anisotropic and very sensitive 
to the scattering geometry when the spiral magnetic or- 
der is the ground state. Our results on the polarization 
dependence of the spectrum suggest that Raman spec- 
troscopy might be very useful to resolve the ambiguity 
in the interpretation of the neutron scattering experi- 



ments and to gain insight into the magnetic structure of 
Cs 2 CuCl 4 . 
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